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f*^ Abstract 

, Using lattice QCD simulation in the quenched approximation, we 

study the re meson, which is 3 Po in the quark model, and compare ex- 

' perimental and other lattice data. The re is the lowest scalar meson with 

i— H , strangeness and constitutes the scalar nonet. The obtained mass is much 

f») ' higher than the recent experimental value, and therefore the re(800) is dif- 

| ficult to consider as a simple two-body constituent-quark structure, and 

, may have another unconventional structure. 

1 Introduction 

ft 

We are in the age of renaissance of hadron spectroscopy, initiated by the an- 
nouncement of the pentaquark baryon pQ , which is followed by the discovery of 
many other possible exotic hadrons with a mass larger than 2 GeV containing 
heavy quarks [2 . These experimental developments prompted the intensive the- 
oretical studies of QCD dynamics with new as well as old ideas on the structure 
and dynamics of the exotic hadrons, such as chiral dynamics [3 , multi-quark 
states with diquark correlations or molecular states and hybrids [2J. 

Such a controversy on the structure of hadrons is also the case for the scalar 
mesons below 1 GeV: the existence of the 1 = and J PC = ++ meson, i.e., 
the cr(400 — 600), has been reconfirmed [H[S] after around twenty years not only 
in 7T7T scattering but also in various decay processes from heavy-quark systems, 
e.g. , D — > tttttt and T(3S) — > Tim [H El [9] . Moreover, the resonance of a 
scalar meson with I = 1/2 is also reported to exist in the K-ir system with a 
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mass m K of about 800 MeV [HI EH EQ • This meson is called the K meson and 
may constitute the nonet scalar state together with the a meson. See FigfT] 



Figure 1: Scalar meson nonet. The a and /o(980) mesons may be ideal mixing 
states of singlet state 4= (uu + dd + ss) and octet state -^=(uu + dd~ 2ss). There 

is , however, experimental evidence that the a meson consist of only uu and dd 
components. Hence, we take the a wave function given in the figure. 

The problem is the nature of these low-lying scalar mesons [12j : they cannot 
be ordinary qq mesons as described in the non-relativistic constituent quark 
model since in such a quark model, the J PC —0 ++ meson is realized in the 3 Pq 
state, which implies that the mass of the a meson must be as high as 1.2 ~ 
1.6 GeV. Thus, the low-lying scalar mesons below 1 GeV have been a source 
of various ideas of exotic structures, as mentioned above: they may be four- 
quark states such as qqqq |13j . or irir or Kir molecules as the recent high- lying 
exotic hadrons can be. These mesons may be collective qq states described as a 
superposition of many atomic qq states |14[ 115) . A mixing with glueball states 
is also possible QH \Tft [HI US] . 

In the previous work [20l [2TI [22l [23] . we have presented a lattice calculation 
for the a meson, by full lattice QCD simulation on the 8 3 x 16 lattice using the 
plaquette action and Wilson fermions: We have shown that the disconnected 
diagram plays an essential role in order to make the a meson mass light. The 
importance of the disconnected diagram suggests that the wave function of the 
a meson may have a significant four-quark, a collective q-q or an even glueball 
component, although the smallness of the lattice requires caution in giving a 
definite conclusion. In contrast to the a meson, the k is a flavor non-singlet 
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state with which a glueball state cannot mix. In previous reports [33[ HI] j we 
reported also a preliminary analysis on the k meson using the dynamical fermion 
for the u(d) quark but using the valence approximation for the s quark, which 
shows that the 1—1/2 scalar meson has a mass as large as about 1.8 GeV and 
cannot be identified with the k meson observed in experiments. 

The lattice volume in the previous investigations was admittedly too small 
to yield a definite conclusion at all, and the lattice cutoff was not appropriately 
chosen to accommodate large masses : m K a > 1, where a is the lattice spacing. 
Hence, we present a simulation with weaker couplings on a larger lattice than 
any other previous simulations although in the quenched level. We perform 
quenched level simulations on the k meson so as to clarify the structure of the 
mysterious scalar meson rather than to reproduce the experimental value of the 
mass; a precise quenched-level simulation should give a rather clear perspective 
on whether the system can fit with the simple constituent-quark model picture 
or not. 

2 Simulation 

We perform a quenched QCD calculation using the Wilson fermions, with the 
plaquette gauge action, on a relatively large lattice (20 3 x 24). 

The values of the hopping parameter for the u/d quark are h u / d = 0.1589, 0.1583 
and 0.1574, while h s = 0.1566 and 0.1557 for the s quark. Using these hop- 
ping parameters except for h s — 0.1557, CP-PACS collaboration performed a 
quenched QCD calculation of the light meson spectrum with a larger lattice 
(32 3 x 56) [55], which we refer to for comparison. The gauge configurations 
are generated by the heat bath algorithm at [3 = 5.9. After 20000 thermal- 
ization iterations, we start to calculate the meson propagators. On every 2000 
configurations, 80 configurations are used for the ensemble average. 

We emply the point-like source and sink for the n + meson 

3 4 
c—1 a—1 

where u{x) and s(x) are the Dirac operators for the u/d and s quarks, and the 
indices c and a denote the color and Dirac-spinor indices, respectively. The 
point source and sink in Eq.([T]) lead a positive spectral function pirn 2 ) in the 
correlation function (k(t)k(Q)) = J dmp(m 2 )eicp(—mt). The result obtained 
here is thus an upper bound of k mass, because our result should include excited 
states. 

First, we check finite lattice volume effects by comparing our results for the 
7r and p masses as well as the mass ratio m^/mp with those of the CP-PACS 
group. The results are summarized in Table [T] Our result for the p meson mass 
is only slightly (< 5 %) larger than the CP-PACS's result. The resulting larger 
value is reasonable because the smaller lattice size gives rise to a mixture of 
higher mass states. We rather emphasize that the deviation between our results 
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and the larger lattice result (CP-PACS) is so small in spite of the large difference 
in the lattice size. 



Table 1: Summary of results for qq type mesons. 



h u /d 


0.1589 


0.1583 


0.1574 


0.1566 


0.1557 


m v 


0.2064(62) 


0.2691(36) 


0.3401(29) 


0.3935(28) 


0.4478(28) 


m p 


0.442(13) 


0.461(06) 


0.496(05) 


0.527(04) 


0.563(03) 


m v /m p 


0.467(21) 


0.584(10) 


0.686(05) 


0.746(03) 


0.796(03) 


m <T v 


1.12(74) 


0.84(23) 


0.886(98) 


0.857(52) 


0.897(35) 






CP-PACS 




25, 


m n 


0.20827(33) 


0.26411(28) 


0.33114(26) 


0.38255(25) 




m p 


0.42391(132) 


0.44514(96) 


0.47862(71) 


0.50900(60) 






0.491(2) 


0.593(1) 


0.692(1) 


0.752(1) 





In Fig. [21 we show m^., ra p and m atl in the lattice unit as a function of the 
inverse hopping parameter l//i u /<j for the u/d quark. The chiral limit (to \ = 0) 
is obtained at h u / d = 0.1598(1) = h crlt (l//i cr j t =6.2581). We find the lattice 
spacing a = 0.1038(33) [fm] in the chiral limit from the value m p a = 0.406(13) 
at this point with the physical p meson mass being used for m p . Note that these 
values are consistent with the CP-PACS's result, h CI n = 0.1598315(68) and a — 
0.1020(8) [fm], within the error bars. 

In Table [T] the mass of the valence a for each hopping parameter is shown; 
the valence a, which is denoted as a v , is defined as the scalar meson described 
solely with the connected propagator. The mass ratio to ct „/to p varies from 2.5 
(hu/d — 0.1589) to 1.6 (h u / d — 0.1557), which is consistent with our previous 
results 23 . In other words, without the disconnected part of the propagator 
the "er" mass becomes heavy. 

The propagators of the K, K* and k mesons are calculated with the same 
configurations using the s-quark hopping parameter, h s = 0.1566 and 0.1557. 
For h s = 0.1557, the effective mass plots of the K* and k mesons are shown 
in Figs. and [31 The masses of the K, K* and k mesons, which are extracted 
from the effective mass plots [26] , are summarized in Tables [2] and [3) Errors 
are estimated by jack-knife method. We find that the effective masses of the K 
and K* mesons have only small errors and are taken to be reliable, while that of 
the k meson suffers from large errors, especially at larger time regions. To avoid 
possible large errors coming from the data at large t, we fit the effective mass of 
the k meson only in the time range 5 < t < 7, 8 where the effective masses are 
almost constant with small errors. Since the effective mass of the K* meson is 
reliable, we show the mass of the k in terms of the ratio to tuk * : Table[4]gives the 
mass ratios vhkI^k- and m K /rriK* at the chiral limit together with m^/niK* 
for h s = 0.1566 and 0.1577. For example, to k /to^« = 0.89(29)/0.4649(69) = 
1.92(61) at h s = 0.1566 in Tabled These calculated mass ratios are shown in 
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Figure 2: mj, m p and m CTl) in the lattice unit as a function of the inverse h u /d- 
The chiral limit is obtained at h cri t = 0.1598(1). 



Fig. All the mass ratios are almost independent of h s . Although the error 
bar for m K /mK* is large, the behavior as a function of h s is reasonable. 

We have searched for the physical value of the s quark hopping parameter 
h s in the following two ways, both of which are found to give similar results: 
1) By tracing a regression line for m^/rriK* (Fig. [5]), we have h s = 0.1563(3) (or 
l/h s = 6.396(13)) for m^/m jK: ,=1019[MeV]/892.0[MeV]=1.143 (input), taken 
from the PDG [3]. This hopping parameter gives the mass ratio m K /mR* — 
1.89(55). 2) We have also determined the hopping parameter so as to reproduce 
the mass ratio m K /m K , = 495.6[MeV]/892[MeV] = 0.5556, with m K = 495.6 
[MeV] being the average value of the Kaon masses given in the PDG [J] . The 
resulting value is found to be h s = 0.1576(2) (or l/h s — 6.3452(80)), which in 
turn gives the mass ratio m K /mK' = 2.00(80). The mass ratios obtained using 
methods 1) and 2) are also presented in Table [4] Both methods give almost 
identical results for the masses of the K, that are about twice that of the K*. 



3 Concluding remarks 

The motivation of our lattice study is to reveal the nature of the scalar meson 
nonet, and the results should be important especially in clarifying how the k 
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Figure 3: Effective mass plots of K* for s quark hopping parameter h s = 0.1557. 



meson with a reported low mass ~ 800 MeV obtained from experiments can be 
compatible with the valence or constituent quark model: the k is a P-wave qq 
bound state in the non-relativistic quark model, and the k meson constitutes a 
nonet together with the a meson and the clq mesons. 

There have not been many lattice studies of k meson. Recently, estimations 
of the k meson have been reported by two groups. Prelovsek et al. [37] have 
presented a rough estimate of the mass of the n as 1.6 GeV, which is obtained 
using the average quark mass of the u and s quarks from the dynamical sim- 
ulations with the degenerate Nf — 2 quarks on a 16 3 x 32 lattice. Mathur et 
al. have studied us meson with the overlap fcrmion in the quenched approxi- 
mation and obtained a mass of the us scalar meson to be 1.41 ± 0.12 GeV [28] , 
The UKQCD Collaboration has studied to some extent the n meson using the 
dynamical Nf=2 sea quarks and a valence strange quark on a 16 3 x 32 lattice 
[29] ; they estimated the k mass as about 1.1 GeV, which is much smaller than 
those in [H [H H7J but still far from the experimental value ~800 MeV. 

In this paper, we have presented the lattice simulation results in the quenched 
approximation for the k meson; the results on the tt, p, K and K* mesons are 
also shown for comparison. 

We have first checked that the masses of the 7r, p, K and K* mesons obtained 
in our simulation are in good agreement with those on a larger lattice (32 3 x 
56) [25] : our results are only within five percent larger than the latter. Our 
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Figure 4: Effective mass plots of the k meson for the s quark hopping parameter 
h. = 0.1557. 



estimated value of the mass of the k is ~ 1.7 GeV, which is larger than twice 
the experimental mass ~ 800 MeV. This result was expected on the basis of 
our experience in calculating the a meson. The relatively heavy mass of the k 
may be at least partly attributed to the absence of the disconnected diagram in 
the k propagator; the k propagator is composed of only a connected diagram. 
While the disconnected diagram was essential for realizing the low-mass a [23] , 
it does not exist for the k; therefore, the mass of the k is not made lighter by the 
disconnected diagram. Indeed, the mass of the valence a v described solely with 
the connected propagator is far larger than the experimental value ~ 500-600 
MeV, as seen in Table [TJ 

Our lattice study and the quark model analysis [30] suggest that the simple 
two-body constituent-quark picture of the k meson does not agree well with the 
experimentally observed k. Note that the quench simulation is a clean theoret- 
ical experiment in which a virtual intermediate like qqqq is highly suppressed 
[13] . Therefore, if its existence with the reported low mass is experimentally 
established, the dynamical quarks may play an essential role for making the k 
mass so lighter or the k may contain an unconventional state such as a qqqq|31| 
or Kir molecular state |32j. which are missing in the calculation here. 

In order to establish this possible scenario, the systematic errors should be 
much reduced in future simulations. Our statistics here is reasonably high (80 
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Table 2: Summary of results for the K, K* and k mesons at h s — 0.1566. 



h u /d 


h 1] . 

"crit 


0.1589 


0.1583 


0.1574 


m K 
m K * 
m K 


0.2829(23) 
0.4649(69) 
0.89(29) 


0.3138(33) 
0.4821(57) 
0.88(23) 


0.3368(30) 
0.4941(49) 
0.81(12) 


0.3677(29) 
0.5117(42) 
0.814(81) 






CP-PACS 






m K 
m K , 




0.30769(28) 
0.46724(84) 


0.32833(26) 
0.47749(74) 





V fecdt = 0.1598(1). 

Table 3: Summary of results for the K, K* and k mesons at h s — 0.1557. 



h u /d 


"crit 


0.1589 


0.1583 


0.1574 


m K 


0.3188(25) 


0.3474(31) 


0.3684(29) 


0.3971(28) 


m K . 


0.4835(61) 


0.5006(52) 


0.5126(44) 


0.5299(37) 


m K 


0.89(21) 


0.88(16) 


0.828(96) 


0.833(72) 



h ciit = 0.1598(1). 



configurations separated by 2000 sweeps) , and the standard meson masses have 
small error bars; see Fig(3] On the contrary, as seen in FigfJ] the effective mass 
of k suffers from large errors for large t, which may be due to a small overlap 
of the physical states. This is not surprising because n is a P-wave meson, 
and expected to be extended. Choosing more adequate extrapolation operators 
and with much higher statistics, we can study the dynamics of hadrons by 
comparing results in the quenched lattice QCD, full lattice QCD and various 
effective theories/models that include the constituent quark models with and 
without the tetra-quark structure, chiral effective theories. 

Acknowledgment T.K. is supported by Grants-in-Aid for Scientific Re- 
search from the Ministry of Education, Culture, Sports, Science and Technol- 



Table 4: Summary of results for the mass ratios uik /mK* an d m K /mK* together 
with m^/rriK* at chiral limit for u/d quarks. 



h s 


0.1566 


0.1557 


0.1563(3) 


0.1576(2) 


l/h s 


6.3857 


6.4226 


6.396(13) 


6.3452(80) 




1.135(10) 


1.164(10) 


1.143^ 




m K /m K , 


0.6086(79) 


0.6593(63) 


0.623(11) 


0.5556 1 ) 


ni K /m K » 


1.92(61) 


1.84(43) 


1.89(55) 


2.00(80) 



^ inputs for calculation of physical value of h s . See the text. 
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